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Abstract: 

We perform the OPE analysis of the contribution of colour-singlet compound states of reggeized 
gluons to a generic hard process in QCD and calculate the spectrum of the corresponding higher 
twist anomalous dimensions in multi-colour limit. These states govern high energy asymptotics 
of the structure functions and their energies define the intercept of the Regge singularities both in 
the Pomeron and the Odderon sectors. We argue that due to nontrivial analytical properties of the 
energy spectrum, the twist expansion does not hold for the gluonic states with the minimal energy 
generating the leading Regge singularities. It is restored however after one takes into account 
the states with larger energies whose contribution to the Regge asymptotics is subleading. 
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1. Introduction 



Scale dependence of hadronic cross-sections is driven by anomalous dimensions calculable in 
perturbative QCD as series in the coupling constant. The classical example is provided by the 
deeply inelastic scattering of a virtual photon 7*(?) off a (polarized) hadron with momentum p^. 
In that case, the operator product expansion (OPE) allows one to expand the moments of the 
structure function F(x, Q 2 ) (with x = Q 2 /2(p ■ q) and Q 2 = —q 2 ) in inverse powers of hard scale 
Q and identify the expansion coefficients as the forward matrix elements of Wilson operators 
.(0) of increasing twist n > 2 and Lorentz spin (j — 1) > 1 

F(j,Q 2 )= / dxxi- 2 F(x,Q 2 ) = J27^ E C nO',«,(Q 2 ))(p|^,i(0)b), (1.1) 

J ° „=2 °* a 

where C^(j, a s ) are the corresponding coefficient functions and the superscript a is introduced to 
enumerate operators of the same twist n. Their total number depends on n and rapidly grows for 
n > 3. The Q 2 — dependence of F(j,Q 2 ) follows from the dependence of the twist— n operators 
on the normalization scale fi 2 = Q 2 . In general, these operators mix under renormalization. 
Diagonalizing the corresponding mixing matrix, one can construct multiplicatively renormalizable 
operators 

Q 2 -^M°n,M\p) = rtU) (pI^(o)Ip) , (I- 2 ) 

where the anomalous dimension has a perturbative expansion = J2'kLi1knU)( a s(Q 2 )/' 7r ) k 

with the expansion coefficients j^nU) having a nontrivial j— dependence. 

For j — > 1, the moments (jl.lj) receive the dominant contribution from the small— x region, 
in which the structure function has a Regge behaviour F(x,Q 2 ) ~ (l/x) a_1 , with a close to 
unity. To find the scale dependence of the structure function F(x, Q 2 ) at small x, one has to 
analytically continue the anomalous dimensions 7^(j) from positive integer j > 2 to "unphysical" 
j ~ 1 and invert the moments (jl.lj) . For twist two this can be done using the well-known DGLAP 
expressions for 7^ =2 (j)- For twist n > 3, the calculation of is much more involved already 
for positive integer j > 2 mainly because the number of operators increases with the twist 
and the size of the corresponding mixing matrices depends on j. As a consequence, analytical 
continuation of anomalous dimensions of the operators of twist n > 3 within the conventional 
OPE approach turns out to be an extremely difficult (if not impossible) task. 

Another approach to finding the asymptotic behaviour of the anomalous dimensions 7^(j) for 
j — > 1 has been proposed in Ref. PQ. It relies on the relation between the twist expansion of the 
moments (jl.lj) and the small— a; behaviour of the structure functions F(x, Q 2 ) obtained within 
the framework of the BFKL approach [2] . Since the anomalous dimensions do not depend on the 
choice of the scattered particles, one can simplify the analysis by considering the deeply inelastic 
scattering of a virtual photon 7*(Q 2 ) off a perturbative "onium" state of the mass p 2 ^ = M 2 , such 
that Aq CD M 2 < Q 2 . In the BFKL approach, one obtains the small— x behaviour in this 
process by resumming a special class of perturbative corrections enhanced by powers of a s ln(l/ x) . 
Among them there are a s \n{l/x) \n(Q 2 /M 2 )— terms responsible for anomalous Q 2 — dependence 
of F(x, Q 2 ). Expanding the resulting expression for the moments F(j, Q 2 ) in powers of 1/Q, one 
can separate the twist— n contribution and calculate the corresponding anomalous dimensions 
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7n(j) f° r J ~~ * I 1 - Still, the actual form of the underlying twist— n operators remains to be found. 
It has been conjectured |H E] that they belong to the class of quasipartonic operators. 

To begin with, we recall that in the BFKL approach, in the generalized leading logarithmic 
approximation (GLLA) j^j, the structure function takes the following form at small— x 

F(x, Q 2 ) = «f ~ 2 Fn{x, Q 2 ) , F N (x, Q 2 ) = J2^M^ BEN{q) ^ (Q)P P q (M) , (1.3) 

N>2 q 

with a s = a s N c /Tr. Here Fn(x,Q 2 ) describes the contribution of color-singlet compound states 
built from N — 2, 3, . . . reggeized gluons, or briefly the N— reggeon states. These states satisfy 
a Schrodinger-like equation for the system of iV gluons interacting on the two-dimensional plane 
of transverse coordinates [BJ. Their spectrum is labelled in (jl.3|) by the set of quantum numbers 
q. For x — > 0, the sum over q in (jl.3j) is dominated by the state with the minimal energy E^iq)- 
The contribution of this state to Fn(x,Q 2 ) has the Regge form F N {x,Q 2 ) ~ (\/ x )-*sE N {q) and 
its Q— dependence is carried by the residue factors (3^,{Q)[3^{M) which measure the overlap of 
the wave function of the state with the wave functions of scattered particles. 

For a s lnx ~ 1 and x — > 0, the sum over N in (jl.3|) is dominated by the N = 2 term and 
gives rise to the BFKL Pomeron [2|. Expanding its contribution to (jl.lj) in powers of 1/Q, one 
finds that the N = 2 reggeon states generate an infinite tower of composite gluonic operators of 
increasing twist n > 2 0j. At n = 2, the anomalous dimensions of twist-two (gluonic) operators, 
7 2 (j), satisfy in the leading logarithmic approximation, a s /(j — 1) = fixed and j — * 1, the master 
equation [I] 

(-^yj = 2^(1) - V(72(i)) - ^(1 - 7 2 (J')) , (1-4) 
where ip(x) = d\nT(x)/dx is the Euler digamma function. Its solution looks like 

7 2 (J) = f^l + 2C(3) (-^yj 4 + 2C(6) (-^y) 6 + 0{al) , (1.5) 

with ((k) the Riemann zeta- function. Eq. (|1.5j) is in agreement with the well-known two-loop 
expression for the twist-two anomalous dimension 7^°- The series (|1.5j) diverges for j — 1 = 0(a s ) 
and its radius of convergence is determined by the right-most Regge singularity, 72 (j) ~ v 'j — 32 
with j 2 — 1 + 4a s In 2 being the intercept of the BFKL Pomeron [2] . 

The contribution of the N > 3 reggeon states to (|1.3|) is suppressed by a power of the coupling 
constant and it serves to unitarize the asymptotic behaviour of the structure function F(x, Q 2 ) at 
very high energy. To perform the OPE analysis of these states, one has to expand the moments 
of Fn(x, Q 2 ) in inverse powers of Q and match the resulting expression for F^(j, Q 2 ) into (jl.lj) . 
Similar to the N = 2 case, the resummed anomalous dimensions of these operators can be deduced 
from the spectrum of the underlying iV— reggeon states. This program has been carried out in 
Refs. P]H| only for the N = 2 states while for N > 3 the main difficulty was a poor understanding 
of the properties of the iV— reggeon states. Later, the anomalous dimensions of twist-four gluonic 
operators have been calculated in the double-logarithmic approximation in Refs. [7J |H]. These 
studies revealed that due to the presence of Regge cuts, the anomalous dimension of N— reggeon 
states has a complicated analytical structure which get simplified however in the multi-color 

x It is known 3], that the OPE expansion breaks down for very small values of x, away from the region 
ln(l/:r) < ln 2 (g 2 /A| CD ). We shall assume that Q 2 is large enough so that this condition is fulfilled. 
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limit. Inspired by this observation, we perform in the present paper the OPE analysis to the 
N > 3 compound states in the limit iV c — > oo. Namely, we identify the leading twist contribution 
of these states to the moments and calculate the corresponding anomalous dimensions for 
j — > 1 in the multi-color limit. Our analysis relies on the remarkable integrability properties of 
the N— reggeon states EI] and on the exact expressions for the energy spectrum of these states 

found in [HI EE!. 

The paper is organized as follows. In Sect. 2 we summarize the properties of the multi-reggeon 
states, perform the OPE expansion of their contribution to the structure function and establish 
the master equation for anomalous dimensions of high twist operators. This equation involves 
analytical continuation of the energy spectrum of the N— reggeon states in multi-colour QCD, 
which is performed in Sect. 3. In Sect. 4 we present our results for the anomalous dimensions. 
Sect. 5 contains concluding remarks. 



2. Reggeon compound states in multi-colour QCD 

Let us perform the OPE expansion of the contribution of the N reggeon compound state to the 
structure function (jl.ljl . Going over to the moment space we find from (jl.3)l 

where the sum runs over the spectrum of the iV— reggeon states. These states satisfy the 
Schrodinger equation which possesses, in the multi-color limit, "hidden" conserved charges q 
and, as a consequence, is completely integrable [HI CHj • 

Due to complete integrability, the spectrum of the N— reggeon states in multi-colour QCD 
is uniquely specified by a complete set of the quantum numbers q = (q 2 , q 2 ■ ■ ■ , qN, Qn) de- 
fined as eigenvalues of the corresponding integrals of motion. The wave function of the state, 
^/ q ({z}; z ), depends on two-dimensional impact parameters of N reggeons, {z} = (zi, . . . , z N ), 
and the impact parameter of the center-of-mass of the state, Zq. Introducing auxiliary complex 
(anti)holomorphic coordinates, z = x + iy and z = x — iy, one finds [3] that the effective QCD 
Hamiltonian for the system of N— reggeons is invariant under conformal, SL(2, C) transforma- 
tions on the plane z — > (az + b)/ (cz + d) and z — > (az + b)/ [cz + d) with ad — be = ad — be = 1. 
As a consequence, its eigenstates ^ q ({z}; z ) belong to irreducible representation of the SL(2, C) 
group: 

N 

V q ({z}; z ) - (cz + d) 2h (cz + d) 2h ]J(cz k + d) 2s (cz k + d) 2s ^ q ({z}; z ) , (2.2) 

k=l 

where (s = 0, s = 1) is the SL(2, C) spin of a single reggeon and the parameters (h, h) define the 
SL(2,C) spin of the iV— reggeon state. For the principal series of the SL(2,C), their possible 
values are parameterized by nonnegative integer and real 

h = — h %v h , h= — h %v h . (2.3) 

By the definition, the wave function ty g ({z}; z ) has to diagonalize the integrals of motion in 
the z— and z— sectors, q = (q2,Q2 ■ ■ ■ ,Qn,Qn)- Also, it has to be a single-valued function on 
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the two-dimensional z— plane and be normalizable with respect to the SL(2, C) scalar product. 
These requirements lead to the quantization conditions for the integrals of motion. 

As was shown in ^2] , the quantized values of q depend on real v% and integer > defined 
in ()2.3|) and on the set of integers £ = (£i, . . . , (.2N-i)- Then, the sum over the iV— reggeon states 
in (j2.1|) looks like Yl q = J2i J^x> & v h Y^ h =o- The eigenstates with different quantum numbers 
q = q(vh,nh,£) are orthogonal with respect to the SX(2,C) invariant scalar product 

N 

/ n rf2 ^(^'^)(^Km;4)r = ^ {2) (^o-4)V, (2.4) 

J k=l 

where 5 qq r = 5{uh — ^'h)5n h n' $w- The energy of the iV— reggeon state, En = EN{vh', n h ■,&)■, is a 
real continuous function of v^. For different integer and £, these functions form an infinite set 
of "trajectories" (see Fig. 1 on the left). 

As we will see in a moment, to calculate the anomalous dimensions we will have to perform 
an analytical continuation of the energy spectrum of the N— reggeon state from "physical" values 
of the conformal spin h, Eq. ()2.3j) . to arbitrary complex h. Since h depends on two parameters, 
rih and one has to decide which of these parameters (if not both) can be made complex. To 
this end, we note that rih and Vh have a different physical meaning. According to (|2.2j) . they 
define the two-dimensional Lorentz spin, h — h = rih, and the scaling dimension, h+h = l + 2ii>h, 
of the iV— reggeon state, respectively. For the wave function ty q ({z}; zo) to be single-valued, rih 
has to be integer while reality condition for follows from the requirement that ty q ({z}; zq) has 
to be normalizable with respect to the scalar product (|2.4j) . Performing analytical continuation, 
we shall require that ty q ({z\] z ) has to be single-valued on the z— plane for arbitrary complex 
h and lift the normalizability condition. This implies that has to be integer and analytical 
continuation goes in u h . The function ^ q ({z}; z ) defined in this way obeys ()2.2jh diagonalizes 
the integrals of motion q but it does not satisfy (|2.4p for arbitrary complex h. 

To expand (|2.1j) in inverse powers of Q, one examines the expression for the residue factors. 
They describe the coupling of the iV— reggeon states to the scattered particles and are given by 

(3 q T {Q) = J d 2 z (%*\V q (z )) , f3 q p {M) = J d 2 z (* q (z )\%) . (2.5) 

Here the scalar product is taken with respect to (|2.4|) while \I/ 7 * ({,?}) and ty p ({z}) are the wave 
functions of the photon and the onium state, respectively, in the impact parameter representation. 
The Zq— integration in ()2.5|) ensures that the total momentum transferred in the t— channel equals 
zero. Neglecting the running of the coupling constant, one finds that the impact factors depend 
on a single scale - the invariant mass of the particle. From dimensional counting, the dependence 
is fixed by the scaling dimension of the iV— reggeon state 

P*.(Q) = C^Q- 1 - 2 ^ , f3*(M) = C q M~ 1+2 ^ , (2.6) 

where C q * and are dimensionless coefficients depending on the charges q = qiyh] nh, £). They 
are different from zero provided that the reggeon state has the same quantum numbers (Bose 
symmetry, C-parity, two-dimensional angular momentum) as hadronic states to which it couples. 
In particular, C q * vanishes for the reggeon states with the Odderon quantum numbers. Our 
subsequent consideration does not rely on the properties of C q * and it is valid for a generic short- 
distance dominated process which receives a nonvanishing contribution from the reggeon states 
both in the Odderon and the Pomeron sectors. 
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Substituting (j2.6j) into ()2.1|) . one gets 



£ n,,>0' 



C«Cl [MY 1+2iUh , s 



Let us examine (|2.7j) in two limits: (i) Q 2 /M 2 = fixed, j — > 1 and (ii) Q 2 /M 2 — > oo, j = fixed. 

In the limit (i), one obtains from (|2.7|) that the rightmost Regge singularity of Fiy(j,Q 2 ) is 
located at — 1 — o7 s min (? .E7v(qr) 



^^'-gi j-j^^ loj ~Q^T (2 - 8) 

Here the dispersion parameter cr^v describes the v^— dependence of the energy E^(q) around its 
minimal value. Eq. (|2.8jl leads to a power rise of the structure function at small— x, F^(x, Q 2 ) ~ 
Q^ 1 (l/x) JJV /(ln 1/x) 1 / 2 , in agreement with the properties of the BFKL Pomeron (N = 2) j2] and 
higher A^ > 3 reggeon states [T2j, 

In the limit (ii), we choose j > and expand the r.h.s. of (|2.7|) in powers of M/Q. For 
M/Q — > 0, the Vh— integration in (|2.7jl can be performed by deforming the integration contour 
into the lower half-plane and picking up the contribution from singularities of the integrand in 
(12. 7|) . The latter could come from singularities of the impact factors C^* , possible branch 
cuts of the energy E^{q) and zeros of the denominator. Let us examine these three possibilities 
one after another. 

It is known 1\ that the impact factors may have poles in due to mixing between gluonic 
and quark operators. Since the mixing does not affect the leading j — > 1 asymptotics of gluonic 
operators, we can safely neglect it. Next, one has to examine the possibility that the energy E N (q) 
has cuts on the complex v^— plane [3 El- in that case, as we will show below, the integration 
around the cut in (|2.7J) provides a nontrivial contribution to F^(j, Q 2 ) which contradicts (jl.lj) . 
For the OPE expansion (jl.lj) to be valid, the contribution of cuts should cancel in the sum 
(12 .7j) over the N— reggeon states. We will argue in Sect. 4 that this is exactly what happens for 
N > 3 reggeon states. Finally, the poles of (j2.7j) originating from zeros of the denominator can 
be defined as solutions to the master equation 

j - l + a s E N {q) = 0, (2.9) 

where q = qiyn] nh,t). At N = 2 this equation matches (jl.4j) for 72 (j) = 1/2 — iv^. For > 3, 
the general solution to (j2.9j) takes the form v h = v h (a s /(j — l);rih,£). As follows from (j2.7jl . 
its contribution to the moments of structure function scales as -F/v(j, Q 2 ) ~ Q~ 2 (M / Q)~ 1+2tVh . 
It matches the OPE expansion, Eqs. (jl.lj) and (jl.2j) . provided that n — 2 / ~f n (j) = 1 + 2iuh, or 
equivalently 

7n(j) = (n - l)/2 - %v h = [ n -{h + h)}/2 , (2.10) 

where h and h are the SX(2,C) spins defined in (12. 3j) . This equation establishes the relation 
between the anomalous dimension of the twist— n operator 7 n (j), Eq. (jl.2j) . and solutions to 
(12 .9j) . Since this operator is built from A^ gluon strength tensors and carries the two-dimensional 
Lorentz spin > 0, its twists satisfies n > N + nh- In addition, writing the anomalous 
dimension as 7 n (j) = 7n°^a s / (j — 1) + 0{a 2 ) 1 one expects that solutions to (j2.9j) should look like 
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ii>h — (n — l)/2 — '-fn^ag/ (j — 1) + C(o!g) with n > iV + n/j. Combining this expression together 
with (j2.9jl one finds 

T (0) 

— r iw9 - (2 - n) 

zz/ fe - (n - l)/2 

Thus, in order for the moments ()2.7|) to admit the OPE expansion (jl.lj) . the energy of the 
iV— reggeon state has to contain an infinite set of poles located along imaginary v^— axis at 
(half)integer points = —i(n — l)/2 with n > N + n^. For iv^ = (n — l)/2 + e the Laurent 
expansion of the energy around the pole at e = can be written as 

E N (q) = — -Co- Cl e + .... (2.12) 

e 

Its substitution into (j2.9j) yields the following expression for the anomalous dimension of the 
twist— n operators, J n (j) — — e ) 



Inti) = c-i 



— + c — +(c -cic_i) — +... 



(2.13) 



with = j — 1. The coefficients depend on the twist of the operator and on the quantum 
numbers specifying the N— reggeon state, = Ck(n,n h ,£). The 0(a s )— term in the r.h.s. of 
(12.13(1 defines the anomalous dimension in the double-logarithmic approximation and 
leads to the following scaling behaviour of the structure function (jl.lj) at small— a; 

F^(x,Q 2 ) ~ ^exp (2 v /c_ 1 ln(l/x)lng 2 ) . (2.14) 

We conclude that in order to calculate the twist— n anomalous dimensions (j2.13j) one has to 
analytically continue the energy of the N— reggeon compound states E N (q) from "physical" values 
of the conformal SL(2, C) spins (12.3(1 to the complex v^— plane and calculate the coefficients of its 
Laurent expansion around the poles located at ivh = (n — l)/2 with n > N + rih, or equivalently 
h — (n + nh)/2 and h — (n — rih)/2. 



3. Analytical continuation 

The N— reggeon states are the eigenstates of the effective QCD Hamiltonian, which is a Hermitian 
operator on the space of functions satisfying (j2.2(l . (j2.3jl and (j2.4(l with Vh real and rih integer. 
The energies of these states, E^{q), are smooth real functions on the real v^— axis and their 
wave functions ty q ({z}] z ) are orthogonal to each other with respect to the scalar product (j2.4jl . 
Performing analytical continuation of the energy spectrum from the real v^— axis to the complex 
Vh— plane, one replaces the normalization condition (j2.4jl by a weaker condition for ty q ({z}; Zq) 
to be a single-valued function on the z— plane. This condition ensures that the two-dimensional 
integrals entering (j2.5jl are well-defined for complex v^- 

From point of view of quantum mechanics, the problem amounts to finding analytical prop- 
erties of the energy E = E(g) as a function of the coupling constant g = v^. In past, this 
fundamental problem has been thoroughly studied in various models [TU EH HZj- It was 
found that, in general, E(g) is a multi-valued function of complex g and the number of its dif- 
ferent branches is equal to the number of the energy levels at real g when the Hamiltonian is 
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hermitian. To study a global behaviour of the function E(g), it proves convenient to glue together 
different "sheets" corresponding to its branches and consider E(g) as a single-valued function on 
the resulting Riemann surface. Depending on the model, this surface may have a rather compli- 
cated structure and consist of a few disconnected parts due to some additional symmetry [To] . 
A remarkable property of the Riemann surface is that it encodes the entire spectrum of the 
model. Namely, knowing the energy of the ground state at real g, one can reconstruct the whole 
spectrum of the model by going around the branching points to other sheets of the Riemann 
surface corresponding to excited energy levels. We shall demonstrate below that upon analytical 
continuation to the complex Vh the energy of the N— reggeon compound states En shares the 
same properties (for iV > 3). 

To begin with, let us consider the well-known expression for the energy of the N = 2 states j2] 

E 2 (v h , n h ) = V + iv^j + $ - iv^J - 2^(1) , (3.1) 

with nh nonnegative integer. E 2 (i'h, n h) is a smooth even function on the real v^— axis. It 
takes its minimal value at Vh — and increases monotonically for Vh — > ±oo. After analytical 
continuation, E 2 (vh,nh) becomes a meromorphic function of Vh- It has an infinite set of poles 
located along imaginary Vh— axis at ivh = ±(n — l)/2 with n > 2 + n^. The leading twist 
contribution, n — 2, corresponds to nh = 0. To obtain the anomalous dimension of twist two, 
Eq. (jl.5|) . one matches the Laurent expansion of E 2 (yh, 0) around the pole at = —i/2 into 
(|2.12jl and applies 1)2.13)1 . The BFKL Pomeron is located on the same complex curve E 2 (vh, 0) at 
Vh = so that E 2 (0, 0) = —4 In 2. For rih > 1 the complex curve E 2 (vh,nh) describes subleading 
Regge singularities and, at the same time, it generates contribution of higher twist. 

For N > 3 reggeon states, analytical continuation of the energy spectrum is more subtle. 
Firstly, the energy of the N— reggeon states does not admit a simple representation like 1)3.1)) 
and, secondly, En = En{v}i', n h, £) depends on the set of 2(iV — 2) integers £, which parameterize 
eigenvalues of the integrals of motion. Thus, performing analytical continuation of En we have 
to deal with an infinite set of complex curves EN(vh',rih,£) labelled by £— integers. For Vh real, 
they define the energy of "physical" N— reggeon states, whose wave functions diagonalize the 
integrals of motion q and satisfy the normalization condition (J2.4j) . Going over to complex Vh, 
one preserves the former condition and relaxes the latter one. 

For Vh real, the integrals of motion in the holomorphic and antiholomorphic sector, and 
q~k, respectively, are conjugated to each other with respect to the SX(2, C) scalar product ()2.4j) 
so that qk = ql (with k = 2, . . . , N). Going over to complex v h , one finds that % 7^ ql and, in 
general, the quantum numbers in the two sectors are independent on each other. Nevertheless, 
the iV— reggeon spectrum contains the states, whose integrals of motion take either real, or 
imaginary values on the real Vh— axis, that is q~k(vh] n>h,£) = ±Q'fc(i / / l ; n>h,£)- For such states, the 
same relation between q^ and % also holds for complex Vh but the charges take complex values. 

Due to complete integrability, the energy of the iV— reggeon state is a function of the integrals 
of motion En = EN{q 2 , q 2 , ■ ■ ■ , <27v, (In) with q^ = qk(vh] nh,£) and similar for %■ For Vh real, the 
energy spectrum is described by infinite set of smooth real functions EN(vh,nh,£) labelled by 
integers and £. The minimal value of En on the real v^— axis determines the position of the 
dominant Regge singularity jn, Eq. (|2.8|) . For Vh complex, one expects that, similar to the N = 2 
case, the complex curve EN(vh] nh, £) has an infinite number of poles located at ivh = ±(n — l)/2 
with n > N + nh- As we will show below, this turns out to be the case but in comparison with 
the N = 2 case one encounters a novel phenomenon. For iV > 3, in addition to the poles, the 
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complex curve E^{vh\ nh, £) contains (an infinite number of) square-root branching points on the 
complex h>h— plane. Therefore, Em is a multi-valued function on the complex v^— plane and its 
different branches are enumerated by integer and £. This property is in agreement with the 
previous findings of Refs. [7J |B|. 

To construct the complex curve E N {v h ; n h , £), we apply the approach developed in Ref. IT2"1. 
Namely, we start with the expression for the energy En obtained there for real Uh and analytically 
continue it to complex following the procedure described above. In this way, one gets 

E N = \ [s(h, q) + e(h, -q) + (e(l - h\ q*))* + - h*, -q*))*] , (3.2) 

where the SL(2, C) spins h and h are given by (J2.3|) and the notation was introduced for the 
quantum numbers in the two sectors, q = {q^} and q = {q^} (with k = 2,. . . ,N). In similar 
manner, -q = {{-l) k j k }, q* = {q* k } and -q* = {(-1)^}, so that E N (-q, -q) = E N (q,q). For 
real v h one has 1 — h* = h and q* = q, so that ()3.2|) produces real values for E N (q,q). The 
function e(h, q) entering ()3.2j) is defined for arbitrary complex h and q as 

e(h, q) = ij e In [e N Q(i + e; h, q)] _ , (3.3) 



e=0 



where Q(u; h,q), the so-called chiral Baxter block, has the following integral representation 

Q(u;h,q)= [ dz z^Qtiz) . (3.4) 
Jo 

Here the function Q\{z) is the solution to the iVth order Fuchsian differential equation 

Qx(z)=0, (3-5) 



N 

\N-k 



(zd z ) N z + (zd z ) N z- 1 - 2(zd z ) N - J2 i k Qk(zd z 

k=2 

with the prescribed asymptotic behaviour at regular singular point z — 1, Qi(z) ~ (1 — z)~ h ~ 1 . 
Since Qi(z) ~ ln^" 1 z for z —> 0, the chiral block Q(u;h,q) defined in ()3.4j) is a meromorphic 
function of u with the iVth order poles located at u = ik for k positive integer. It is convenient 
to normalize Q\(z) in (j3.4j) in such a way that the residue at the iVth order pole u = i equals 
unity. Together with (j3.3J) this leads to 

To evaluate e(h,q), one has to solve the differential equation ()3.5|) and match the resulting 
expression for Q(u;h,q), Eq. ()3.4j) . into ()3.6|) . This allows one to determine ()3.2j) for arbitrary 
complex q and q. 

The quantization conditions for the integrals of motion q = (q 2 , q 2 , ■ ■ ■ , qN, ^Af) follow from 
the requirement for their common eigenfunctions ^ q ({z}; Zq) to be single-valued functions on the 
z— plane. These conditions take the form ^2] 

Q(i + e; h, q)Q{i - e; h, -q) -Q(i + e;l- h, q)Q(i -e;l-h, -q) = 0(e°) . (3.7) 

Substituting ()3.6|) into this relation, one finds that the l.h.s. scales as l/e 27V_1 for e — > 0. Require- 
ment for the poles at e = to vanish leads to the overdetermined system of 2N — 1 equations 
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for the 2(N — 2) charges qt and q~k (k — 3, . . . ,N). The two remaining charges are given by 
q2 = —h(h — 1) and q~2 = —h(h — 1). Solving this system, one obtains the quantized values of the 
integrals of motion and verifies that the obtained expressions satisfy three consistency conditions. 

Applying Eqs. ()3.2|) - ()3.7|) one can calculate the energy of the N— reggeon state for arbitrary 
complex v h . For v h real, they lead to the results for the energy E N {y h \n h ,t) obtained before 
in Ref. ^2]. For Vh complex, they define analytical continuation of the function EN{vh',nh,£)- 
At iV = 2 the differential equation (|3.5j) can be solved in terms of Legendre functions and the 
resulting expression for the energy ()3.2j) coincides with the known expression (|3.1j) . For N > 3 
the analysis is more involved since ()3.5|) cannot be solved exactly and one has to rely on a power- 
series solutions to (|3.5|) described at length in Ref. [TS]. In the next section, we summarize the 
results for iV = 3, 4, 5, 6 reggeon states. 

4. Spectral surface 

At iV = 2, the energy spectrum E2(vh,nh) is described on the complex z/^— plane by the family 
of meromorphic functions (J3.1)) labelled by a nonnegative integer rih- For iV > 3, an analytical 
expression for the energy is not available but the value of En can be calculated for arbitrary 
complex Vh using Eqs. ()3.2j) - (j3.7|) as follows. Let us consider an arbitrary contour on the complex 
Vh— plane that starts on the real axis and terminates at some complex v^. For given real Vh, 
the energy E N (v h ; n h ,£) takes a discrete set of real "physical" values labelled by integer n h and 
£ = (Ei, . . . , £2n-a) (see the left panel in Fig.[T} [HI [12]. Calculating the energy along this contour 
point-by-point from Eqs. ()3.2j) - ()3.7j) . one analytically continues the function EN(vh',nh,£) to 
complex Vh. 

This procedure allows us to determine the global properties of the complex curve E^iyh] rih, £) 
on the complex Vh— plane. If En were a single-valued function of complex Vh, it would resume 
its original value after going around arbitrary closed contour on the v^— plane. This is the case 
at N = 2, Eq. (J3.1)) . whereas for N > 3 one finds that E^iyh] n h, is a multi-valued function of 
Vh- Namely, E^ has (an infinite number of) branching points on the complex plane, Vh = Vbr,k, 
such that En changes its value after encircling these points. We found however that En resumes 
its value if one encircles the branching point twice. This implies that En has square-root cuts 



where E N defines the energy on the upper and lower edges of the cut, respectively. Aside of the 
branch cuts, En has an infinite number of poles. In agreement with our expectations, Eq. (j2.11|) . 
they are located along imaginary axis at ivh = ±(n — l)/2 with n > N + n^. 

Thus, for N > 3 the energy E^ivh] nh, £) is a meromorphic function on the complex Vh— plane 
with the square- root cuts running between the branching points zv,fc- Another way to represent 
the family of multivalued functions EN{vh'i n h-,£) is to sew together their branches along the 
square- root cuts, Eq. (|4.1jl . and define E N as a single- valued meromorphic function on the result- 
ing Riemann surface. Following [TJ], we shall call it the spectral surface. Its topology depends 
on the number of reggeons N. Different sheets of this surface can be enumerated by integer rih 
and £. 

Let us consider a particular (n h ,£)—th sheet of the spectral surface lying over the v h — plane 
and suppose that it is sewed with (n' h ,£')—th sheet along the square-root cut that starts at 
the branching point iv- For real Vh, the value of function En on these two sheets defines the 




(4.1) 
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"physical" energy of two iV— reggeon states, E^iyh', rih, £) and Ex{yh\ n' h ,£'), respectively. If one 
analytically continues these functions along the contour that starts at real Vh and terminates at 
the branching point, = iv> then at the vicinity of the branching point the functions behave 
as (j4.1|) leading to 

E N (u hl ;n h ,£) = E N (u hr ; n' h , i ') . (4.2) 

Thus, in a complete analogy with quantum mechanics ^3], the branching cuts (|4.1J) arise due 
to collision of the energy levels at some complex Uh = away from the real v^— axis. Since 
the wave functions of the two states have to coincide at the branching point ^3], they have to 



have the same two-dimensional Lorentz spin, 



n 



hi 



and possess the same quantum numbers 



(quasimomentum, C-parity, Bose symmetry, etc) leading to additional selection rules for £ and 
£'. This implies that the spectral surface can not be simply connected and it should consist 
of (infinite number of) disconnected components enumerated by nonnegative integer rih and the 
quantum numbers just mentioned. 



Ert 




Figure 1: The energy spectrum of the N = 3 reggeon states E^{vh]Uh,£) for rih = and £ = 
(0,^2)) with £2 = 2, 4, . . . , 14 from the bottom to the top (on the left). Analytical continuation 
of the energy along the imaginary v^— axis (on the right). The branching points are indicated by 
open circles. The lines connecting the branching points represent Re E3. 



To illustrate our findings, let us consider the N = 3 reggeon states with the Lorentz spin 
rih = 0. Their spectrum is specified by two integers £ = (£1,^2)- Assigning the values of 
these integers to different energy levels, we follow the convention adopted in Ref. The 
energy E 3 (v h ; 0, £) is a smooth even function on the real v h — axis (see Fig. [T] on the left). It 
approaches the minimal value min^-E^ = .24717 at = along the "trajectory" with l\ — 
and £ 2 = 2. The charges qziyh) an d q~z{vh) take pure imaginary values along this trajectory, so 
that ^3 + q3 = 0. After analytical continuation, the same relation holds for arbitrary complex 
including the branching point ()4.2j) . Therefore, the energy level with l\ = and ^2 = 2 could 
only collide with the levels for which 53 + ^3 = 0. As was shown in ^2]; the latter condition is 
satisfied for l\ — and £ 2 = even positive. Examples of such energy levels with £ 2 = 2,4,..., 14 
are shown in Fig. ^ They define a connected component of the N = 3 spectral surface. 

To elucidate analytical properties of the energy, we used Eqs. (I3.2|) - (|3.7j) to analytically 
continue eight functions E3(vh]0,£ = (0, £2)) with £ 2 = 2, 4,..., 14. These functions define eight 
different branches of the complex curve E$(vh) which can be represented as a two-dimensional 
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surface in C 2 = M 4 space with the (uh, E$) —coordinates. The slices of this surface along Im Vh = 
and Rez/fc = hyperplanes are shown in Fig. Q on the left and the right panels, respectively. We 
observe that, firstly, the energy has poles at iz/ pole = 3/2, 5/2, 7/2, 9/2, 11/2 and, secondly, 
different energy levels collide at iv hr = 2.70, 3.29, 4.73, 5.34, 5.74, 5.94. At the vicinity of the 
branching point the energies of two colliding levels behave as E s , Eq. (|4.1|) . As expected, the 
poles of E3 are located along the imaginary z/^— axis at iuh = (n — l)/2 with n = 4, 6,.... This is 
not the case, however, for the branching points. It turns out that some of the branching points 
are located away from the imaginary v^— axis and, therefore, they can not be seen in Fig. ^ For 
instance, the branches with £2 = 2 and £2 — 4 collide at iu^ = 1.723 + i .248. In other words, the 
"ground state" branch, £ 2 = 2, is sewed with the £2 = 4 branch which in its turn is sewed with 
the £2 = 6 branch and so on. Thus, going around the contour on the complex Uh— plane which 
starts at some real Vh, encircles the branching points and returns to the starting point one can 
reconstruct the energy spectrum of the N = 3 reggeon states with £\ = and £2 = even positive. 

The fact that the branching points take complex irrational values implies that the con- 
tribution of the corresponding square-root cuts to the integral ()2.7|) scales as l/[Q 1+2tUbt In 3 / 2 Q] 
and, therefore, it breaks the OPE expansion (jl.lj) . Although such corrections are present in ()2.7|) 
for given nh and £, they cancel against each other in the sum over all states. To see this let 
us examine two terms in the sum (|2.7j) corresponding to the energy levels colliding at Uh = Vbr- 
Each of them is associated with a particular sheet of the spectral surface. Denoting the energy 
on the upper (+) and the lower (— ) edges of the square-root cut on these two sheets as Ef{y h ) 
and E^(uh) one finds that Ef(v h ) = E^{u h ). As a consequence, the contour integrals around 
the same cut on the two sheets differ by a sign and their sum equals zero. Thus, the contribution 
of the cuts to ()2.7)1 cancels completely in the sum over all N = 3 reggeon states belonging to 
the same spectral surface. At the same time, if one had retained in (j2.7j) only contribution of 
the state with the minimal energy E^ t) - the one defining the leading Regge singularity (j2.8j) . 
the OPE expansion would have been broken. To restore the OPE one has to keep in (|2.7|) the 
contribution of subleading Regge singularities. 

As another important example, we consider analytical continuation of the N = 3 reggeon 
state with the Lorentz spin = 1, the so-called descendant state [T2j. This state has q^Vh) = 
on the real v^— axis and its energy coincides with the energy of the N = 2 state, E 3 ^ = E^ivn'-, 1), 
Eq. (|3.ip . Both relations survive the analytical continuation and hold true for arbitrary complex 
Vh. Thus, in distinction with the previous case, the energy of this state E 3 ^(i/h) is a single-valued, 
meromorphic function on the complex v^— plane. At Vh = it takes the value E% = £2(0; 1) = 
which defines the "physical" ground state energy for the system of N = 3 reggeons ^H]- Going 
over to the imaginary z/^— axis, one finds that E 3 has poles in the lower half-plane at iv^ = 1, 2, ... 

We would like to stress that the N = 3 states with = and n h = 1 described above 
define two different components of the unifying N = 3 spectral surface. The reason why we 
have selected them is that they contain two special states E^{vh = 0;n^ = 1) = and E^iv^ = 
0; nh = 0) = .24717, which are currently considered as solutions for the Odderon state in QCD 
with the intercept jo = 1 — otgE^ [T9*l ITS] . As was shown in Sect. 2, the poles of the energy in 
the lower half-plane Im Uh < 0, induce the contribution to the OPE expansion (jl.lj) of the twist 
n = h + h = l + 2ivh. The minimal twist n min corresponds to the pole closest to the origin. 
As follows from our analysis, n m \ n = 3 for the descendant state with rth = 1 2 and n m ; n = 4 

2 We disagree on this point with Ref. [5], in which it was claimed that n m i n = 4 for the Odderon state with the 
intercept 1. 
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for the state with rih = 0. Thus, the leading twist contribution of the two above mentioned 
Odderon states with the intercepts jo = 1 and jo < 1 scales at large Q 2 as ~ 1/Q 3 and ~ 1/Q 4 , 
respectively. 



5. Anomalous dimensions 

To calculate the anomalous dimensions (j2.13j) . one has to work out the Laurent expansion of 
the energy around its poles, Eq. (|2.12|) . We recall that for the N— reggeon state the twist— n 
contribution comes from poles located at iv^ = (n — l)/2 with n > N + n^, or equivalently 
h — (n + rif l )/2 and h = (n — rift)/2. 

Let us first consider the N = 3 descendant state, nu = 1. For complex its energy, E 3i d, 
coincides with E^iyhi 1) defined in (jH.lj) . The closest to the origin pole is located at ivh = 1 and 
the corresponding SL(2, C) spins (|2.3|) are given by h = 1 + iv h = 2 and h = iv h = 1 so that the 
twist equals n = h + h = 3. Using (|3.ip . one finds the expansion of the energy E 3<d = E 3 d (h + e) 
for e — > as 

£ 3 , d (2 + e) = - + 1 - e - (2C(3) - 1) e 2 + . . . (5.1) 
e 

Applying Eqs. ()2.12|) and (|2.13|) . one finds from (|5.1|) the twist-3 anomalous dimension corre- 
sponding to the Odderon state with the intercept j© = 1 as (u = j — 1) 



(N=3) 

7s 



(j) = ^-p) 2 + (2C(3) + l)p) 4 + ..., (5.2) 



where the subscript and the superscript indicate the twist and the number of reggeons entering 
the state, respectively. 

Let us now consider the N = 3 states with rih = 0. Their analytical properties are encoded 
in the spectral surface shown in Fig. [TJ Applying Eqs. (|3.2|) - (l3.7j) . we calculate the energy at 
the vicinity of the poles at h = 1/2 + iv^ = 2, 3, 4, 5 and obtain the following expressions for 
E 3 = E 3 (h + e) 

E 3 (2 + e) = £-! + !- I e + 1.7021 e 2 + . . . 

£ 3 (3 + e) = 2e" 1 + — -1.6172 e + 0.719 e 2 + ... 

8 

(4 + e) = e~ l + 0.6806 e - 1.966 e 2 + . . . 

12 

^ b) (4 + e) = 2e- 1 + ^-3.2187e + 3.430 e 2 + ... 

^ a) (5 + e) = 2e^ + — - 2.0687 e + 1.047 e 2 + . . . 

d v ' 48 

^ b) (5 + e) = 2 e" 1 + — - 2.4225 e + 0.247 e 2 + . . . (5.3) 

Here ellipses denote 0(e 3 ) terms and the additional superscript was introduced to distinguish 
between different branches. We recall that all states in ()5.3j) have the same Lorentz spin = 0, 
so that h — h — 1/2 + u>k and the corresponding twist equals n = 2h > 4. Applying Eqs. ()2.12|) 
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and (|2.13|) . one finds from the first relation in (|5.3)l that the leading, twist-4 anomalous dimension 
corresponding to the Odderon state with the intercept jo < 1 is given by 

7r 3, ( . ) = |_l ( ^_l ( | )3 _ L077l( ^ + ... (, 4) 

The remaining relations in (J5.3)) lead to similar expressions for higher-twist anomalous dimen- 
sions. To save space we do not present them here. The following comments are in order. 

The first few terms of the Laurent expansion in ()5.3|) can be calculated exactly. The reason for 
this is that solutions to the differential equation ()3.5|) can be expanded in powers of e and the first 
few terms can be obtained in a closed form. The coefficient in front of 1/e in ()5.3|) equals either 
1 (only for h even), or 2. In the latter case, writing the energy as E^(h + e) = 2e _1 + 7(/i) + 0(e) 
one finds that finite O(e )— terms have the following remarkable property. It turns out that 
7(4) = 15/4 and 7(5) = 53/12 coincide with the energy of the Heisenberg SL(2, R) magnet model 
of spin s = 1. Most importantly, the energy spectrum of this model determines the anomalous 
dimensions of local composite three-quark (baryonic) operators of helicity— 3/2 [201 E] For 
such operators, integer (h — 3) counts the number of covariant derivatives and 7(^1) defines their 
anomalous dimension. Similar relation holds for higher N between the energy of the N— reggeon 
states (SL(2, C) Heisenberg spin magnet) and anomalous dimensions of N— particle operators 
(SX(2,R) Heisenberg spin magnet). Its origin will be discussed in a forthcoming publication. 

Going over to the N > 4 case, we shall concentrate on the N— reggeon states providing the 
contribution to the OPE expansion (jl.lj) of the minimal twist n m i n . As before, we shall denote 
the corresponding anomalous dimension as 7n2 in (j). 

For N = even we find that the minimal twist equals the number of reggeons involved n min = N. 
It corresponds to the pole of the complex curve En located at iuh = (N — l)/2 and rih = 0, or 
equivalently h = h = N/2. This pole is situated on the same spectral surface as the "physical" 
ground iV— reggeon state, that is the state with the minimal energy for real v^. Its energy, 
E^ = mm Uhnh £ E N (is h ; n h} £), defines the intercept of the N— reggeon states in the Pomeron 
sector [II1II2]: e[ st) = -.67416 and E^ v) = -.39458. The expansion of the energy E N (h + e) 
around e = looks like 

2 1 

£4(2 + e) = - + 1 - -e - 1.2021 e 2 + . . . 

£ 6 (3 + e) = - + ^--^ e - 0.238e 2 + ... (5.5) 
e 2 16 

Applying Eqs. ()2.12j) and (|2.13jh one obtains the following expressions for the leading, twist— N 
anomalous dimension of the N— reggeon states 

7 r 4) (j) = 2^ + 2(^) 2 -13.6168(^) 4 + ... 

T f =6) (j) = 4^ + 6^) 2 + 2^) 3 - 33.23 ^) 4 + ... (5.6) 

UJ \U J \UJ J \UJ J 

with uj = j — 1. We observe that higher order corrections to (J5.6j) are large indicating that 
the series have a finite radius of convergence and its value decreases with N. This is in a 
qualitative agreement with the fact that the intercept of the N— reggeon states scales at large 
iV as \]n — 1| ~ 1/N [12J. Writing the 0(a s )— correction to (j5.6j) as 72n(^) = ^n«s/^ (with 
n = 2,3), one finds that e n verifies the condition e n < 2{n — 1) established in [T3]. In addition, 
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72n(<^) < n72(co , /n) which means that the anomalous dimension of iV— reggeon states is subleading 
in the multi-color limit as compared with the anomalous dimension of N/2 BFKL Pomerons (ZlEj- 
For iV = odd, similar to the N = 3 case, we consider separately the sectors with the Lorentz 
spin rih = and = 1. In the first case, the minimal twist equals n m j n = (N + 1) and it 
corresponds to the pole of the energy at iv^ = N/2 and = 0, or equivalently h = h = (N+l)/2. 
For instance, at iV = 5 the expansion of the energy E^{h + t) and the anomalous dimension look 
like 

3 7 101 „ , , 

£ 5 3 + e = - + --_ e -0.1136e 2 + ... . 5.7 
e o 21b 

ir%) = 3^ + !(5i) 2 -I(5i) 3 - 13.032(5i)V... 

u 2 \ uj J 8 V w / \uj J 

The pole (|5.7j) is located on the same spectral surface as the physical ground state in this sector, 
E^ gr \n h = 0) = .12751. This state has the quantum numbers of the Odderon. It has the intercept 
1 — a s E^ which is smaller then 1 and its leading twist equals 6. 

In the second case, for iV = odd and = 1, the minimal twist is smaller n m i n = N. It 
corresponds to the pole of the energy En at ivh = (N — l)/2, or equivalently h = (N + l)/2 and 
h = (N — l)/2. At N = 5 the expansion of the energy E^{h + e) looks like 

£ 5jd (3 + e) = 3 + — + 1.36180 - 0.4349 e - 0.315 e 2 + . . . (5.8) 

7 f = 5 >(j) = 3 -±^^ + 3.56524f^) 2 + 1.8743f^) 3 - 11.21 9 f^) 4 + . . . 

As in the iV = 3 case, this reggeon state is descendant [12] , that is its energy equals the energy of 
the N = 4 state with rih — 1. The pole ()5.8|) belongs to the same spectral surface as the N = 5 
ground state, E ( 5 SI) {nu = 1) = 0, which has the quantum numbers of the Odderon and intercept 
equal to unity. The second relation in ()5.8|) defines the leading, twist— 5 anomalous dimension of 
this state. 



6. Conclusions 

In this paper, we performed the OPE analysis of the contribution of iV— reggeon compound 
states to the moments of the structure function F(j, Q 2 ) for j — > 1 and calculated the anomalous 
dimensions of the leading twist contribution in multi-color QCD. To this end, we analytically 
continued the energy of the iV— reggeon states E^ from the "physical" region of parameters 
(real v^— axis) to complex Vh and established the relation between the anomalous dimensions 
for j — > 1 and the Laurent expansion of En around its poles. At N = 2 the energy is a 
meromorphic function on the complex v^— plane [3] while for iV > 3 analytical properties of the 
energy E^ are changed dramatically. Namely, we found that, in agreement with previous findings 
[3IH], the energy En is a multi- valued function on the complex v^— plane. The reason for this 
is that different energy levels in the spectrum of the N — reggeon states collide after analytical 
continuation to complex and, as a consequence, their energies develop square-root cuts. Due 
to nonvanishing contribution of the cuts, each energy level (the sheet of the spectral surface) 
breaks the twist expansion of F(j, Q 2 ), but it is restored in the sum over all A— reggeon states. 
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To summarize our main results, we found that the leading contribution of the iV— reggeon 
states to the moments (jl.lj) has the twist N. For even and odd iV it comes from the reggeon states 
with the two-dimensional angular momentum = and = 1, respectively. For N = 3, 4, 5, 6 
the corresponding anomalous dimensions j£ are given by Eqs. ()5.2j) . ()5.6|) . ()5.8|) . The explicit 
form of the underlying twist— N operators need to be found. We demonstrated that two solutions 
for the Odderon states with the intercepts jo = 1 and jo < 1, recently discussed in the literature, 
have the twist 3 and 4, respectively. Their anomalous dimensions are defined in Eqs. (|5.2|) and 
(|5.4jl . Notice that one-loop corrections to the twist— N anomalous dimensions, Eqs. ()5.6j) - (j5.8j) . 
have a rather simple form. It would be interesting to compare these expressions with those 
obtained by explicit calculation of the twist— iV anomalous dimensions in the double- logarithmic 
approximation (see e.g. [T3]). 

Finally, let us comment on the relation of our results to those obtained in Ref. 0. In that 
paper, the anomalous dimension of the N = 3 and N = 4 reggeon states have been calculated 
in the multi-color limit and the expressions obtained there differ from Eqs. (|5.2|) and ()5.6|) . The 
reason for disagreement is the following. The approach proposed in [5] is based on the assumption 
that for rih ^ the energy of the N— reggeon states as a function of conformal spins h = (l+7i/i)/2 
and h = (1 — n/J/2 can be obtained by analytical continuation in h from the energy in the sector 
with the Lorentz spin nh = 0. As follows from our analysis, this assumption is erroneous. For 
N > 3 the energy is a multi-valued function of conformal spin (see Fig. ^) and it can only 
be analytically continued over the spectral surface described in Sect. 4. The states with different 
spins rih belong to different disconnected components of this surface so that their energies are not 
related to each other by analytical continuation. If one ignored nontrivial analytical properties 
of the energy and applied the approach of Ref. p. , one would generate spurious states which do 
not belong to the physical spectrum of the reggeon states. 
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